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ABSTRACT 


necent  developments  In  the  theory  of  the  sine-Gordon  equation 


are  used  to  analyze  the  appearance  of  a "displaced  linear  branch"  in 


the  volt-ampere  characteristic  of  a large  Josephson  Junction.  Internally 


the  junction  is  divided  into  a flux  annihilation  domain  near  the  center 


and  flux  flow  domains  near  the  edges.  The  displaced  linear  branch  is  expert 


mental  evidence  for  the  existence  of  flux  flow  domains.  In  the  flux  annihilation 


domain,  an  ac  voltage  component  is  induced  by  the  continuous  formation 


and  decay  of  sine-Gordon  breathers.  Whitham's  nonlinear  WKB  method 


is  used  to  analyze  the  flux  flow  domains  while  the  inverse  scattering 


transform  method  is  used  as  a qualitative  guide  to  the  study  of  breather 


formation  in  the  flux  annihilation  domain 
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MAGNETIC  FLUX  ANNIHILATION  IN  A LARGE  JOSEPHSON  JUNCTION 

Alwyn  C.  Scott 


I.  INTRODUCTION 

The  physical  phenomenon  to  be  discussed  in  this  paper  is  the 
"displaced  linear  branch"  in  the  volt-ampere  characteristics  of  Josephson 
tunnel  Junctions  having  dimensions  large  compared  with  the  Josephson 
penetration  length  0.1  mm).  This  effect  was  reported  several  years 
ago  [1-3]  and  an  explanation  based  on  a continuous  flow  of  magnetic 
flux  quanta  (fluxons)  was  suggested.  Since  that  time,  several  striking 
advances  have  been  made  in  the  understanding  of  the  analytic  behavior 
exhibited  by  solutions  of  the  underlying  nonlinear  partial  differential 
equation,  the  sine -Gordon  equation 

*xx  ' *tt  = Sln  * * (L1) 

Particularly  important  have  been  i)  the  development  of  an  inverse 
scattering  transform  method  (ISTM)  for  (1.1)  [ 3 - 6] , ii)  the  application 
of  Whitham's  method  (WM),  of  slow  perturbation  of  a single  phase 
(periodic  traveling  wave)  solution,  to  a dissipative  version  of  (1.1)  [7,  8] 

4*xj{  ‘ ♦tt  " a\  = Sin  *»  (1,  2) 

ill)  Ablowitz's  discussions  of  the  extension  of  Whitham's  method  to 
multiple  phase  (multiply  periodic)  solutions  [9  - 12] , and  iv)  demonstrations 
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that  appropriate  multiply  periodic  solutions  assumed  by  Ablowitz  do  in 
fact  exist  for  (1.1)  [13  - 15] . The  present  discussion  should  be  con- 
sidered in  context  with  previous  efforts  to  apply  the  fruits  of  soliton 
theory  to  the  problems  of  fluxon  dynamics  in  large  Josephson  Junctions 
[16,  17] . Such  studies  are  of  interest  not  only  as  an  aid  to  understanding 
the  volt-ampere  characteristics  of  Josephson  junctions  [1  - 3,  18,  19] , but 
also  because  the  soliton  properties  of  fluxons  may  make  them  useful  as 
carriers  of  information  [ 20  - 25] . 

The  displaced  linear  branch  was  first  observed  on  1 mm.  x 1 mm. 
crossed  strip  Josephson  junctions  as  volt-ampere  characteristics  of  the 

sort  displayed  in  Figure  1.  Above  a critical  current  (I  ) for  zero  voltage 

c 

(-  1/2  ampere),  the  voltage  (V)  was  found  to  Increase  linearly  with  the 
difference  between  current  and  critical  current.  Thus 


V = K(I-I  ) 
c 


(1.3) 


where  the  constant  of  proportionality  (K)  is  much  less  sensitive  to 
temperature  than  the  normal  electron  component  of  tunneling  current.  It 
was  suggested  [1]  that  the  input  power  (VI)  implied  by  (1.  3)  is  carried 
by  moving  fluxons  from  the  edges  of  the  Junction  toward  the  center 
where  it  is  given  up  in  events  of  fluxon- a ntifluxon  annihilation. 
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1.35  *K 


(b) 


(c  ) 


Figure  1.  Voit-ampere  characteristics  of  typical  experimental  junctions 
at  various  temperatures,  (a),  (b),  and  (c)  are  Sn-Ox-Sn,  vertical:  wOO 
mA/div,  horizontal:  2 mV/div  (d),  (e),  and  (f)  Pb-Ox-Pb,  vertical  500 
mA/div,  horizontal:  2mV/div. 


3- 


mm 


Although  the  two  dimensional  effects  which  appear  in  the  square 
junction  geometry  of  Figure  1 are  of  both  theoretical  and  practical  interest, 
the  present  discussion  will  be  confined  to  the  one  dimensional  geometry 
indicated  in  Figure  2.  In  this  case  a mechanical  analog  of  (1.  2)  [ 26]  can 
be  used  to  demonstrate  the  existence  of  flux  flow  domains  near  the  edges 
and  a flux  annihilation  domain  near  the  center  as  shown  in  Figure  3. 

The  paper  is  organized  in  the  following  way.  The  linearized  version 
of  (1.  2)  which  is  obtained  by  neglecting  the  sin  <j>  term  is  discussed  in 
Section  II  in  order  to  provide  a perspective  for  considering  the  nonlinear 
problem.  Then  using  results  obtained  by  Pelinovskii  and  Shavratzkii  [8] 
and  in  [ 17]  the  flux  flow  domain  is  described  in  Section  III.  The 
nonlinear  solution  in  the  flux  flow  domain  is  based  upon  the  single 
phase  (periodic  traveling  wave)  solution  described  in  detail  by 
Whitham  [7]  . The  flux  annihilation  domain  considered  in  Section  IV 
involves  the  propagation  of  fluxons  in  both  directions  and  requires 
therefore  a double  phase  (doubly  periodic)  solution  [9  - 12]  . In  the 
flux  annihilation  domain,  fluxons  and  antifluxons  combine  to  form 
"breather  solitons"  which  give  up  their  energy  to  the  dissipation 
represented  by  the  -ad>t  term  in  (1.  2).  To  understand  the  annihilation 
domain,  it  is  important  to  appreciate  the  nature  of  breather  solitons; 
thus  two  appendices  are  included  which  discuss  (A)  the  structure  of  a 


I *********  • 


Figure  2.  One  dimensional  geometry  used  for  the  theoretical  analysis  (not 
to  scale),  (a)  perspective  view,  (b)  section  showing  voltage  v across 
insulating  barrier  and  current  i parallel  to  the  barrier. 


breather,  and  (B)  the  use  of  the  inverse  scattering  transform  method 
(ISTM)  to  study  the  generation  of  breathers  through  the  interaction 
of  double  phase  components.  A third  appendix  (C)  gives  a summary 
account  of  single  phase  and  double  phase  solutions  for  (1.1).  The 
material  in  these  appendices  should  be  appreciated  in  order  to  follow 
the  discussion  of  the  main  text.  Finally  we  conclude  with  a review  of 
the  results  which  appear  to  be  useful  in  interpreting  currently  available 
experimental  data  and  some  suggestions  for  future  experiments. 

The  unnormalized  form  of  (1.  2)  is 

*xx  - 1 g*TT  - Q*  *T  = V sin(2ir«/®0)  U*  4) 

where 

X is  distance  in  the  laboratory  coordinate  system, 

T is  time  in  the  laboratory  coordinate  system, 

9 = V(X,T)  is  the  voltage  across  the  barrier, 

S>  = -I(X,  T)  is  the  superconducting  surface  current  flowing 
parallel  to  the  barrier, 

g is  the  shunt  conductance  through  the  barrier  carried  by 
(Giaever  type  [27] ) normal  electron  tunneling  per 
unit  length, 

i is  the  inductance  presented  to  the  superconducting  surface 
current  per  unit  length , 

c is  the  shunt  capapacitance  of  the  barrier  per  unit  length, 
and  JQ  is  the  maximum  (Josephson  type  [28])  superconducting  electron 
tunneling  through  the  barrier  per  unit  length. 
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Thus 


•2 


(X,  T)  * / V(X,  T’)dT' 


(1.5) 


is  the  magnetic  flux  which  has  passed  point  X at  time  T.  It  is 
convenient  to  measure  this  in  units  of  the  magnetic  flux  quantum  or 
fluxon 


4>0  = h/2e 


15 


(1.6) 


= 2.067843  X 10  volt-seconds. 


If  we  define  a Josephson  length 


V 


<2 


1/2 


(1.7) 


and  note  that  the  characteristic  velocities  of  (1.4)  are  ± uQ  where 


“„  = l<c]-1/2, 


(1.*) 


then  (1.4)  reduces  to  (1.2)  with  the  normalizations 


x = XA 


I 


t = u0TAj 


(1.9) 

(1.10) 


<j>  = 2 w*/4>0 


« = gl  UQ  . 


(1.11) 

(1.12) 
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(1. 15a,  b) 


In  this  paper,  we  will  be  primarily  concerned  with  an  analysis 
of  the  normalized  equation  (1.  2).  A more  detailed  discussion  of  the 
approximations  involved  in  obtaining  (1.4)  can  be  found  in  reference  [17] 
together  with  some  tables  of  experimentally  measured  values  of  the 


parameters. 


II.  THE  UNEAR  APPROXIMATION 


We  begin  with  a brief  discussion  of  the  equation 


♦xx  • *tt  '*V°  <2-‘» 

which  is  a linearization  of  (1.  2)  by  neglecting  the  sine  term.  From  (1.4) 
it  can  be  seen  that  this  is  physically  equivalent  to  assuming  the  coef- 
ficient of  Joseph  son  tunneling  current  (jQ)  to  be  small.  Equation  (2. 1) 
can  also  be  written  as  the  first  order  pde's 


v + i = -ay 
t x 

(2.  2a, b) 

i + v = 0 
t x 

where  v and  i are  normalized  shunt  voltage  and  series  current  as 
defined  in  (1.13)  through  (1.15).  The  experimental  situation  indicated 
in  Figure  2b  can  then  be  represented  as  in  Figure  4a  where  the  dimension 
b = BAj.  Physical  symmetry  and  the  boundary  conditions  require  the 
current  (i)  to  be  zero  at  the  center  in  Figure  4a.  Thus  we  need  consider 
only  the  half  Junction  shown  in  Figure  4b  with  the  boundary  conditions 


v(0,  t)  = vQ 
i(b/2,  t)  = 0. 


(2.  3a, b) 


If  we  assume  vQ  = const,  and  seek  steady  state  solutions,  (2.2)  implies 
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j'"  >#. 


(2.4a,b) 


v(x)  = vQ 

<*bv„  _ 

1(x)  = “T~  (1 ' T1 


Now  return  to  (1.15)  and  write 


4>  = vt  - ix 


(2.  5) 


which  can  be  viewed  as  a traveling  wave  with  phase  velocity 


u = v/i  . 


(2.6) 


From  (2.  4)  we  see  that 

~ < 1 =o  u > 1 for  0 < x < b/  2 


and 


(2.7a) 


— > 1 =>  u < 1 for  0 < x < l/o 

u > 1 for  l/a  < x < b/  2 . 


(2.7b) 


Thus  for  sufficiently  small  a,  the  phase  velocity  of  a single  traveling 
wave  (2.  5)  will  exceed  unity  for  almost  all  x . 

We  can,  if  we  wish,  resolve  <(>  into  the  sum  of  two  traveling 
waves,  one  carrying  flux  to  the  right  and  the  other  carrying  flux  to  the 
left.  Then 
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(2.8) 


4*  = (Vjt  - ijX)  + (v^t  + i2x)  . 


where  Vj  + = v^  and  ij  + i^  = i(x).  This  resolution  is  not  unique, 

but  in  anticipation  of  results  to  be  obtained  below  for  the  nonlinear  problem 
we  suggest  the  following  for  the  case  ab/2  > 1. 

For  0 < x < l/o  set  = 0 and  i2  = 0 and  take  Vj  and  ij 
to  be  the  values  given  in  (2.4).  For  l/a  < x < b/2  take 


i 


1 


b-l/a-x 

(b7TI7? 


i 


2 


_v  , h ) 

2 " 2 'b/ 2-1/ a 


(2. 9a, b) 


Thus  for  0 < x < l/a  there  is  only  a single  traveling  wave  component 
and  it  has  a phase  velocity  less  than  unity.  For  l/a  < x < b/2  there  are 
two  components  and  the  phase  velocity  of  each  has  a magnitude  equal  to 
unity. 
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HI.  THE  FLUX  FLOW  DOMAIN 

In  this  section  we  address  ourselves  to  the  nonlinear  equation 
(1.2).  Consider  first  a steady  state  solution.  Then  ^ = sin  <t> 
which  is  readily  integrated  to 

- 4>Z  = E - cos  <(>  *3,1) 

2 9x 

where  E is  a constant  of  integration.  The  conditions  4>  - 0 and 

a - 0 at  large  x requires  E = 1 whereupon 
^x 

a = ± 2 sin  <)>/  2 . 

^x 

From  (1.15)  i has  maximum  magnitude  of  ic  = 2 when  4>(°>  = * ‘ Thus 
the  critical  current  (into  both  ends) 


!/2 


r®  i 

„ oo 
Uni  . 


(3.2) 


is  established  below  which  the  voltage  (v  = 4>t>  is  zero  and  above 
which  flux  must  flow.  This  effect  has  been  carefully  investigated  [29]  • 
Suppose  now  that  1(0)  > ic  (I  > Ic>  so  a steady  state  solution 
for  (1.  2)  does  not  exist.  As  has  been  shown  in  reference  [ 17] , it  is 
convenient  to  turn  to  Whitham's  method  of  averaging  for  which  it  is 
supposed  that  the  solution  appears  locally  as  a traveling  wave  of  the  form 


4>  * 4»(© ) 


0 = wt  - kx 


where  w and  k are  allowed  to  be  slowly  varying  functions  of  x 


and  t.  Then  4.(6)  is  defined  locally  as  the  elliptic  integral 

e JJ7J  f , (3.4 

where  we  assume  E > 1.  Dynamic  equations  for  the  slow  variation  of 
k,  w and  E with  x and  t are  [7,  8,  17] 


[ wF'(E)F(E)] . + [ kF' (E)F(E)]  = -au>F'(E)F(E) 

l X 

k + w = 0 . 
t x 


(3.  5a, b) 


The  functions  F'(E)  and  F(E)  are  complete  elliptic  integrals  of  the 
first  and  second  kind  respectively  and  defined  as 


r'(EI*i  /"  75<rhos  *) 

. 2it  

F(E)  = r — f "J  2(E  - cos  4>)  d4- 
2 w J 


(3. 6a, b) 


Note  that  in  the  large  amplitude  linear  limit.  E •*  00 , 


F'  - (2E) 


-1/2 


F - (2E) 


(3.7a,b) 


It  is  convenient  to  require  that  <fi  increase  by  2rr  when  0 does, 
whereupon  k,  « and  E are  related  by  the  nonlinear  dispersion  equation 
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k2  - W2  = [ F'(E)] 


(3.8) 


In  the  large  amplitude  linear  limit  (E  » 1),  the  product 

F'(E)F(E)  - 1 


as  is  shown  in  Figure  5.  Furthermore  (3.4)  and  (3.8)  imply 

4>  — 6 

so 


v = 4>t  et  = w 


(3.9a,b) 


and  (3.  5)  approach  the  linear  equations  (2.  2)  which  were  considered  in 
the  previous  section.  Figure  5 indicates  that  this  large  amplitude  linear 
limit  is  approached  within  a few  percent  for  E>  2. 

Suppose  now  that  the  current  (i)  is  increased  to  the  maximum 
zero  voltage  (v  = 0)  level.  As  i goes  from  0 to  0+(  the  steady 
state  current,  i(x),  will  Jump  from  the  solution  of  (3.1)  (with  E = 1) 
to  a solution  of  (3.8)  with  w = 0 or 

k = k = [ F'(E  )]_1  (3.10) 

c c 


where  E is  the  critical  value  of  the  integration  constant  E for  a 
c 

traveling  wave  solution  carrying  the  critical  current.  Since  k^  c ic  = 2, 

(3.10)  is  satisfied  for  F’(E  ) * 0.  5 or 

c 

E « 2. 237. 
c 

For  positive  values  of  voltage  (w  > 0),  this  periodic  solution  will 
move. 

To  determine  the  relation  between  terminal  current,  i(0),  and 
voltage,  w(0),  whon  the  flux  is  moving,  we  must  eliminate  E from  the 
nonlinear  dispersion  equation  (3.8).  This  can  be  done  by  noting  [7,  17] 
that  traveling  wave  velocity  (or  nonlinear  phase  velocity) 

u = w/k  (3.  U) 

and  E are  related  as 

=/  [--^r- 

0 l-(u')  E 1 ' ' J 

c 

where  the  lower  limits  have  been  chosen  to  satisfy  the  condition 

u = 0 when  E = E . 

c 

1/2 

Since  E * 2,  the  high  amplitude  linear  limit  [ F(E)  « (2E) 
c 

and  F*(E)  * (2E)  -i/a, 

is  appropriate  for  evaluating  the  right  hand  side 

of  (3.12).  Then 
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which  upon  substitution  into  (3.8)  gives 


u>(0)  = k(0)  - siZZ 

o 


(3.14a) 


This  is  the  equation  for  the  displaced  linear  branch.  Taking  account  of 
(3.9)  and  the  normalizations  in  (1.13)  and  (1.14),  it  can  be  written 


V 


i & «- 


y 


(3.14b) 


Thus  the  constant  slope  (K)  which  appears  in  (1.  3)  is  equal  to  half  the 
characteristic  impedance  of  the  linear  operator  in  (1.  2)  [ 1] , and  the 
weak  temperature  dependence  of  this  slope  (see  Figure  1)  is  readily 
appreciated. 

Now  let  us  consider  how  far  such  a steady  state  traveling  wave 
solution  can  be  extended  into  the  Junction.  Uie  assumption  of  steady 
state  in  (3.  5b)  implies 


cu  = const  (3.15a) 

and  from  (3.5a),  k(x)  must  satisfy 

[ kFF'l  = -a-wFF'  . (3.15b) 

1 Jx 

Pelinovskii  and  Shavratzkii  [ 8]  have  shown  that  there  is  a definite 
limit  in  x beyond  which  a single  phase  (periodic  traveling  wave)  solution 


of  (3.15)  does  not  exist.  To  see  this  consider  the  quantity  kFF', 
Using  (3.8)  it  is  seen  to  be 

kFF'  = fVi  + w2(F')2 


(3.16) 


which,  for  fixed  co,  is  a function  of  E with  a minimum  value  (E  = Emin) 
satisfying  the  condition 


“2«r'<Em1n)'2+F(En,ln»F"(E— »tl=0- 


min 


(3.17) 


The  character  of  kFF'  as  a function  of  E and  u is  displayed  in 
Figure  6 and  Emin  vs.  u is  plotted  in  Figure  7.  When  kFF'  is  equal 
to  its  minimum  value,  (3.15b)  can  no  longer  be  satisfied  because  the  right 
hand  side  requires  that  kFF'  continue  to  decrease.  It  is  only  for  the 
range  of  x for  which  E > that  the  steady  state  single  phase 

solution  can  exist. 

The  value  of  x (say  x^)  at  which  kFF’  has  fallen  to  its 
minimum  value  is  an  upper  estimate  for  the  boundary  of  the  flux  flow 
region.  Equation  (3.15)  can  be  integrated  with  the  boundary  conditions 

E = E(0)  at  x = 0 


E = E . at  x = xrt 
min  0 


to  obtain 


2,rn2 


1 + u(F') 


dE  . (3.18) 


where  from  (3.13)  and  (3.14) 


E(0)  = (W2Ec  + Ec 


(3.19) 


1/  2 

In  the  large  amplitude  linear  limit  (3.18)  implies  that  xQ  - (2Ec)  /aw 
as  «-*■<»,  and  xQ  -*■  .607/ aw  as  w -►  0. 

We  are  now  prepared  to  calculate  the  power  flowing  across  the 
boundaries  of  the  flux  flow  region.  In  general  P(x)  = <v(x)i(x)>  * wk(x) 
where  the  approximation  improves  for  large  E as  indicated  in  (3.9). 
Taking  the  boundaries  of  the  flux  flow  region  as  x = 0 and  x = xQ 
(see  Figure  8),  the  power  flowing  in  at  x = 0 is 

P(0)  = wk(0) 


= W2  + oa 1 2E 

c 


where  (3.14a)  has  been  used  to  evaluate  k(0).  In  a similar  way,  the 
power  flowing  out  at  x = x^  is 


P(xQ)  = wk(xQ) 


where  (3.8)  has  been  used  to  evaluate  k(xQ). 


Figure  8,  Steady  state  components  of  voltage  (w)  and  current  (k)  in  the  flux 
flow  domains.  In  the  flux  annihilation  domain,  an  approximation  to  the  double 
phase  solution  is  sketched. 


If  the  assumption  of  steady  state  in  the  flux  flow  domain  is 
correct,  the  difference  between  these  incoming  and  outgoing  powers  must 
equal  the  d.c.  dissipation.  Thus  we  expect 


P(0)  - P(x0)  = au1 xQ 


(3.  22) 


whereupon  a second  estimate  for  xQ  is 


. 1 
x„  = — 
0 a 


1 + 


/ 2E 
v c 


2E 


- /I  + 


min 


(3.  23) 


We  do  not  expect  precise  agreement  between  the  values  of  xQ  calculated 
from  (3.18)  and  (3.  23)  because  the  power  flow  estimates  in  (3.  20a)  and 
(3.  21a)  are  not  exact  for  smaller  values  of  E. 

However,  as  Figure  9 shows,  the  agreement  is  rather  good. 
Equation  (3.23)  gives  values  for  xQ  which  are  7-10%  higher  than  those 
calculated  from  (3.18)  when  w lies  between  1 and  2.  Equation  (3.23) 
is,  of  course,  much  more  convenient  than  (3.18).  For  w > 10,  the 
asymptotic  expression 


(3.24) 


is  a reasonable  approximation. 


IV.  THE  FLUX  ANNIHILATION  DOMAIN 


Let  us  suppose  that  we  have  a large  Joseph  son  junction  which 
exhibits  a displaced  linear  branch  (Figure  1)  in  the  volt-ampere  char- 
acteristics as  is  implied  by  (3.14).  Assume  further  that  a,  b,  w and 

E are  adjusted  such  that  the  relation 

c 

orb  2 r~ 

— w =W(W+V2EC)  (4.1) 

holds.  This  relation  says  that  the  d.c.  power  into  the  Junction  (RHS) 
equals  the  d.c.  power  dissipated  on  the  Junction  (LHS).  If  w is 
increased  above  the  value  which  satisfies  (4.1),  the  LHS  will  be  greater 
than  the  RHS  and  the  displaced  linear  branch  will  not  supply  power  to 
the  junction  at  the  rate  it  would  be  dissipated.  Thus  the  value  of  w 
which  satisfies  (4.1)  is  the  "break  voltage"  u>  (see  Figure  10)  above 
which  the  displaced  linear  branch  will  not  continue.  Figure  1 shows 
that  such  a break  voltage  is  a characteristic  feature  of  the  displaced 
linear  branch. 

For  Wg  » 1,  (4.1)  implies  that  a necessary  condition  for 
observing  a break  in  the  displaced  linear  branch  is 

orb 

~~  1 (4.2a) 

or,  in  unnormalized  units, 
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(4.  2b) 


Now  consider  the  difference  between  the  power  flowing  nto 
the  flux  annihilation  domain,  P(x0>,  anc*  the  d.c.  power  being 
dissipated  in  this  domain 


pdC  ■ ““2<!  - v 


(4.3) 


Thus,  from  (3.21)  and  (3.23) 


P(xft)  - P , = u>2  + u / 2E  - 

0 dc  v c 


2 b 

arw  — . 


(4.4) 


From  (4.1)  this  difference  is  zero  when  u>  = w . For  w < w (4.  4) 

B B 9 

indicates  that  the  d.c.  power  flow  into  the  flux  annihilation  domain  is 
greater  than  the  d.c.  power  dissipation  in  that  domain.  Thus  for  u>  < w 

B 

it  is  not  possible  to  have  a steady  state  (w  = const)  in  the  flux 
annihilation  domain.  There  must  be  an  additional  a.c.  component  of 
voltage  to  dissipate  the  power  difference  given  by  (4.  4).  The  source 
for  this  a.c.  voltage  is  the  continuous  creation  of  breather  solitons  near 
the  center  of  the  Junction. 

The  dynamic  activity  in  the  flux  annihilation  domain  is  consider- 
ably more  complex  than  in  the  flux  flow  domains.  To  get  some  qualitative 
appreciation  for  this  activity  we  turn  to  the  ideas  developed  in  Section  II 
and  in  the  appendices.  Consider  first  the  problem  of  representing  the 


■ 


d.c.  components  or  the  time  averages  of  v and  i.  Following 


Appendix  C we  assume  a double  phase  solution 

4>(x,  t)  = *(0^  e2) 


with 


0j  s w^t  - k^x 


e„  = w t - k*x 
2 2 2 


Then,  motivated  by  (2.8),  we  try  the  simple  form 


®(er  e2)  * eA  + e2  . 


Equations  (C- 14)  both  become 


,k>  tkz'x‘  -“(“i  + "21 


and  (C-6c,  d)  imply 


w + « = w (a  constant). 

The  boundary  condition  kj(Xg)  + k2(xQ)  = anc* 

kl  + k2  = k<x0>  ' ‘ xo>  ' 
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(4.  5) 

(4.6a,b) 

(4.7) 

(4.8a) 

(4.8b) 

(4.  8a)  implies 

(4.9) 


The  symmetry  condition  discussed  in  connection  with  Figure  4 requires 

k](^)  + k^)  = °»  but  from  *4*9) 


kl(l)  + k2(I)  = k(X0)  ‘ aw(2  ' V 


(4.10) 


= u + / 2E  .^u. 
vc  2 


This  is  equal  to  zero  only  at  the  break  voltage  w = w . For  w < w , 

D D 

the  simple  double  phase  function  assumed  in  (4.7)  cannot  correctly 
represent  the  d.c.  components.  Some  more  versatile  expression 
(probably  using  Riemann  theta  functions)  will  be  required  [13  - 15]  . 

In  the  special  case  « = however,  we  can  use  (4.  5)  - (4.7)  where 

k<y  x ■ *<> 

Vx>  = k,V  - “,x-xo’ + ~T  'virr1 


k(x0)  x • xo 

k2(xi — • 


For  w - w « w (4.11)  should  be  approximately  correct. 

D D 


Furthermore  from  (3.  24)  and  (4.  2a) 


so  most  of  the  Junction  will  be  occupied  by  the  flux  annihilation  domain. 


Also  the  rate  at  which  kinks  enter  this  domain  at  x = xQ  will  be 
w/2it.  Since  each  breather  forms  by  coupling  a kink  entering  from  the 
left  and  an  antikink  from  the  right,  the  steady  state  rate  of  breather 
formation  must  equal  w/2rr. 

Where  in  the  flux  annihilation  domain  will  these  breathers  form  ? 

To  answer  this  question  we  can  use  the  results  of  Appendix  B as  a qualitative 

guide.  Just  at  the  center  of  the  Junction,  <j>  =0  and  <*,  = w.  Thus, 

x t 

from  (B-2Z),  we  expect  the  number  of  breathers  formed  per  unit  time  and 
distance  to  be  aw/2ir.  Since  the  rate  of  rotation  falls  from  a maximum  at 
the  center  to  zero  at  the  edges  of  the  flux  annihilation  domain,  the  net 
rate  of  breather  formation  should  equal  abw/4iT  . Then  the  condition 
(4.  2a)  implies  a breather  formation  rate  ~ w/2-rr. 

Thus  we  have  a fairly  self  consistent  picture  of  breather  formation 
for  the  case  « " wg  « WB*  exPect  breathers  to  remain  near  the 
center  of  the  Junction  because  they  are  formed  near  the  center  of  the 
Junction  as  stationary  breathers. 

The  case  w « is  much  less  well  understood  for  several  reasons. 


i) 


The  qualitative  nature  of  the  d.c.  solution  is  not  clear. 


il)  From  (4.4)  the  a.c.  voltage  is  no  longer  a small  fraction 
of  the  d.c.  voltage. 

ill)  For  1,  x0  ~ b/2  and  the  flux  annihilation  domain 
becomes  very  narrow.  In  this  case,  the  a.c.  voltage  is  probably  not 
confined  to  the  flux  annihilation  domain. 

Numerical  studies  of  the  case  w « w would  be  helpful  in 
understanding  these  uncertain  aspects.  As  far  as  the  external  observation 
of  a displaced  linear  branch  is  concerned,  however,  all  that  is  required 
is  for  a single  phase  (traveling  wave)  solution  to  be  established  near  the 
edges  of  the  Junction. 
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V.  CONCLUSIONS 

Here  an  attempt  will  be  made  to  draw  together  the  somewhat 
tattered  strands  of  the  above  discussion  by  indicating  which  results  may 
have  experimental  significance  and  by  suggesting  some  directions  for 
future  research. 

1)  Critical  current. 

From  (3.2),  the  ciritical  current  (I  ) above  which  flux  will  begin 

c 

1/2 

to  move  is  4[4,0J0/2iri  ] . Since  it  is  practical  to  compute  JQ  and  l 

from  Independent  measurements  on  a small  Junction  [17]  , the  ideal  value 

for  I can  be  compared  with  that  actually  observed  on  a large  Junction, 
c 

Such  a comparison  is  important  because  irregularities  in  the  barfier 

thickness  (d  in  Figure  2)  may  cause  flux  pinning  and  increase  I above 

c 

the  value  given  by  (3.  2). 

2)  Slope  of  the  displaced  linear  branch. 

Equation  (3.14b)  indicates  that  the  slope  of  the  displaced  linear 
branch  should  equal  2\lc/l  where  the  factor  of  two  accounts  for  current 
flow  into  both  ends  of  a long  Junction  (see  Figure  2b).  This  result  con- 
firms the  suggestion  made  in  [ 1]  . Since  c and  I can  be  independently 
measured  on  a small  Junction,  this  striking  feature  of  the  above  theory 
for  the  displaced  linear  branch  is  always  subject  to  direct  experimental 
check. 


-34- 


In  (4.  2)  it  was  noted  that  the  condition  for  observing  a break  in 


the  displaced  linear  branch  is  ab  * 2 or,  in  unnormalized  units, 
gB  * 2\c/f  . If  ab  » 2,  the  Josephson  current  becomes  a small 
portion  of  the  total  current,  and  the  situation  reduces  to  the  linear 
problem  studied  in  Section  II.  If,  on  the  other  hand,  ab  « 2,  the 
Junction  becomes  a resonant  electromagnetic  structure,  and  the  step 
structure  observed  by  Fiske  [ 30-33]  will  dominate  the  volt-ampere 
characteristic.  Thus  the  condition 

— ~ 1 ( 5.  la , b) 

or 

gB  ~ Z^c/l 

also  determines  the  range  of  experimental  parameters  for  which  the 
displaced  linear  branch  should  be  unambiguously  observed. 

4)  Numerical  studies. 

It  was  assumed,  in  the  discussions  of  Sections  III  and  IV,  that  the 
a.c.  voltage  associated  with  the  decay  of  breather  solitons  is  confined 
to  the  flux  annihilation  domain  and  does  not  "leak  out"  into  the  flux  flow 
domain.  This  assumption  seems  plausible  because  stationary  breathers 
are  generated  at  the  center  of  the  Junction  and  it  is  confirmed  in  a 


-35- 


rough  way  by  observations  on  the  mechanical  model  as  shown  in  Figure  3. 

But  it  would  be  comforting  to  have  a direct  numerical  test  of  this  picture. 

An  appropriate  problem  would  be  integration  of  (1.  2)  with  end  conditions 

w(0,  t)  = w(b,  t)  = 0 for  t < 0 

= v for  t > 0 

and  with  ab  = 2,  a = 0.1,  0.01,  0.001,  and  v = 1,  10,  100.  An  asymptotic 
form  of  the  dynamic  activity  should  be  sought  as  t 00 . 

5)  Experiments  with  light  sensitive  junctions. 

Barone  et  al.  [ 34]  have  recently  demonstrated  that  properly 
prepared  Pb  - CdS  - In  Junctions  can  be  changed  from  " small"  (b  < 1) 
to  "large"  (b  > 1)  by  ordinary  optical  illumination.  This  change  occurs 
because  the  Josephson  tunneling  is  increased  and  therefore  kj  is  decreased 
under  illumination.  An  increase  in  the  Josephson  tunneling  should  correspond 
to  an  increase  in  the  normal  electron  (Giaever)  tunneling  and  therefore  an 
increase  in  a.  Thus  ab  should  be  an  even  more  sensitive  function  of 
illumination  than  b.  Indeed,  from  (5.16),  ab  is  directly  proportional  to 
the  Giaver  tunneling  density. 

Thus  it  may  be  possible  to  experimentally  observe  a transition  from 
"Fiske  steps"  on  a "small"  Junction  to  flux  flow  on  a "large"  Junction  as 
the  illumination  level  is  increased. 
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APPENDIX  A:  STRUCTURE  OF  A BREATHER 


1.  Classical  structure 

Consider  the  sine-Gordon  equation  in  characteristic  coordinates 


4>£T  = sin  4,  (A-l 

which  is  equivalent  to  (1.1)  under  the  independent  variable  transformation 


% = \ (x-t) 

T - \ (x+t)  . 


(A-2a,b) 


If  is  a solution  of  (A-l),  it  is  easily  shown  that  <J>n(£,T) 

which  satisfies 


'Wi’e 3 4iY 

, * -*  , (A’3o-b| 

....  . i . , n vn-l. 

V*n-l  T = ' 7 3ln(— l » 

is  also  a solution  [ 35] . This  is  called  a Backlund  transform  which 
generates  a new  solution  from  a known  (old)  solution;  by  repeated 
application  a hierarchy  of  solutions  can  be  obtained.  If  i^  is  real, 
an  additional  kink  or  soliton  will  be  introduced  with  the  asymptotic 
form 

<J,j  = 4 tan_1[  exp(2iy^  - iV^)] 
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(A-4) 


and  velocity  in  the  laboratory  (x,  t)  frame  [17] 


4v  +1 


UJ  = 4Vj-l 


As  George  Lamb  has  suggested  [ 36]  , a convenient  way  to  generate  a 
breather  solution  to  (1.1)  is  to  perform  two  Backlund  transforms  with 
complex  y's  and 


Yl  = 


(A-6) 


Such  a solution,  <j>  , is  given  by  [16,  35  , 36] 


* tan  "4 


V*2 


(A- 7) 


where  cfcj  and  4>2  have  the  Guder mannlan  form  Indicated  in  (A-4).  Thus 
in  the  laboratory  (x,  t)  frame  a breather  sollton  has  the  form 


4,  =4 tan  < 


.if/  *\  exp[  l(y - ~ )x-i(y+  ^ )t]  -exp[  -i(y* - )x+l(/+ )t] 


\w+\  / l+exp[i(v-y*  - ^ + ^7‘)x+  K-V+V* -4\  + ^)tl 


4y  4y 


(A-8) 


Upon  examination  of  (A-8),  the  following  points  will  be  noted. 


i)  As  v -*  a pure  imaginary  #,  -*  2*  and  the  breather 

disassociates  into  a kink -antikink  pair. 

ii)  As  v areal#,  <t>B  — 0 and  is  proportional  to 

2 2 

sin(kx-wt)  where  k = y « 1/4  y,  u = y + 1/4  y and  w = 1 + k . 


ih 


iii)  The  envelope  of  the  breather  is  determined  by  the  denominator 


in  (A-8).  The  envelope  velocity 


u = 
e 


. J_  * 1 

V + 4v 


1 


* 4.  1 

V + 


(A-9) 


4y  1 4y’' 

is  positive  for  |y  | > 1/2  and  negative  for  |y  | < 1/2.  Thus  the  condition 
for  a stationary  breather  is 


M = 4 


(A  - 10) 


In  terms  of  an  angle  parameter  (p)  which  relaxes  from  — -*•  0 as  two 
kinds  decay  into  a stationary  breather  of  zero  energy,  (see  Figure  11) 
a stationary  breather  takes  the  form 


4>b  = 4 tan  1 {tan  p sin[  (cos  P)t]  sech[  (sin  P)x]  } 


(A  - 11 ) 


The  rest  energy  of  a breather  can  be  calculated  as  the  total 
energy  of  a stationary  breather  from  the  Hamiltonian  density 


12  12 

H = “ <t>x  + ~ 4>t  + 1 - c°s  4> 


(A -12) 


corresponding  to  (1.1).  At  t = 0,  (A-ll)  reduces  to  <j>g(x,0)  = 0;  so 


1 2 


H = — $8  t(x,0)  and  the  total  energy 

oo 

UB  = / Hdx 


-00 
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i 


w 


t 


or 


UD  = 16  sin  P . 

D 


(A  -13) 


When  p reaches  ir/2,  the  breather  disassociates  into  two  kirks  of 
the  form  (A-4),  so  the  energy  of  a kink  is 


Uk=8 


(A-14) 


2.  Quantum  structure 

Fadeev  [37,  38]  and  Dashen,  Hasslacher  and  Neveu  [39]  have 
carried  through  an  exact  quantization  of  the  breather  soliton  described 
by  (A-ll).  The  discrete  energy  levels  were  found  to  be 

n 


U = 2U.  sin 
n k 


w. 


(n  + 1/2) 


(A-13) 


where  is  the  unnormalized  energy  of  a kink  (fluxon),  is  Planck's 
constant  divided  by  2ir  and 


“I  = U0AJ 


(A-14) 


is  the  Josephson  frequency  which  was  used  in  (1.10)  to  normalize  the  v 
time  scale.  There  is  correspondence  between  the  classical  frequency 
(cos  P)  in  (A-ll)  and  the  frequency  for  emission  or  absorption  in  (A-13) 
If 


k 


(A-15) 
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Thus  the  allowed  quantum  states  for  a stationary  breather  are  evenly 
spaced  on  the  circle  ly  I = 1/2  (see  Figure  11).  Quantum  effects  should 
become  important  in  predicting  the  dynamics  of  a breather  if  the  number 
of  these  states  is  not  large  compared  with  unity.  To  evaluate  this 
situation  we  note  that  in  unnormalized  units 


\ * 8cVnxj  • ; Voxj 


(A-16) 


and  turn  to  the  parameters  measured  on  the  Josephson  transmission  lines 
discussed  in  reference  [17] . From  (A-14)  the  total  number  cf  levels  is 


(A-17) 


Table  I.  Estimates  of  Quantum  Structure  for  Breathers 


X. 

J 

WJ 

N53C 

2.63  X 10'4 

N25L 

1. 27  X 10'3 

meter 

6.69  x 1010 

\10 

1.81  X 10 

rad/sec. 

Jo 

1.9 

.097 

amp/  meter 

Jo 

1. 22 

.062 

amp/CM^ 

uk 

-19 

13.  1 X 10  7 

-19 

3.24  X 10 

Joules 

N(=4Uk/n^t>J) 

2.  37  X 105 

2.  16  X 105 

— 

From  the  last  entry  in  Table  I it  is  clear  that  the  number  of 
quantum  levels  is  so  large  that  classical  dynamics  should  be  sufficiently 


I 

- ( 

accurate  for  the  description  of  these  particular  losephson  transmission 

2 

lines.  As  for  other  structures  it  should  be  noted  that  N <*  JQX-j  or 

N<*W  (A-18) 

where  W is  the  width  of  the  transmission  line  (see  Figure  2a).  If 
W were  decreased  by  two  orders  of  magnitude  (from  64  microns  to 
0.64  micron),  quantum  corrections  would  still  be  of  minor  importance. 
This  is  in  accord  with  a previous  estimate  of  the  need  for  quantum 
corrections  in  a point  contact  Josephson  Junction  [ 40]  . 


APPENDIX  B:  A COHERENT  BREATHER  GAS  FROM  ROTATING 
INITIAL  CONDITIONS 


Here  we  use  the  inverse  scattering  tran^Wir.  method  for  the  sine- 
Gordon  equation  (1.1),  which  was  discovered  byAblowitz,  Kamp,  Newell 
and  Segur  [ 4]  , to  see  how  a breather  gas  can  be  generated.  In  laboratory 
coordinates  the  equations  for  evolution  of  the  scattering  variable 
4»  = (4»i»  ^2)  are 

+1,*  = ■ i • 77 cos  ♦>  h * i 1 ir sln  * - i (vv)4,2 

(B-la,bl 

*2,x  " i [ * S‘"  *+  2 W1*!  + i lv  • 7;  cos  *>+2 


and 


1, 

2,t 


= 2 (v  + cos  ^ ~ Sin  4>  + | (4>x-4>t)H2  . 

= i ^ sln  * - 1 (Wl4,l  ‘ 2 (V  + ~ti  COS  ^2  • 


(B- 2a, b) 


These  equations  are  closely  related  to  the  Backlund  transform  equations 

(A- 3)  [ 5,  41  ] for  which  we  noted  that  a purely  imaginary  y introduced 

* 

a kink  or  soliton  and  a pair  (\,  ) introduced  a breather  into  the 

total  solution. 

The  application  of  such  scattering  equations  to  compute  the 
evolution  of  nonlinear  waves  from  specified  initial  conditions  has  been 
discussed  in  detail  [ 5,  17] . Here  I will  merely  sketch  a particular 
calculations  of  interest  to  our  study  of  the  flux  annihilation  domain  in 

Section  IV. 
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Refer  to  Figure  12  and  consider  the  following  initial  conditions 
on  <t>(x,  t)  which  is  to  be  a solution  of  (1.1) 

4>(x,  0)  = 0 

<j>t(x,  0)  = 0 for  lx  | >p  (B-3a,b) 

= V for  |x  | < p . 

A scattering  problem  is  defined  as  indicated  where  4^  and  4^  are  the 
asymptotic  Incident  and  reflected  waves  in  Regional)  and  4^  is  the 
asymptotic  transmitted  wave  in  Region(T)  . Since  4<  obeys  (B-l,  2)  in 
all  regions,  the  asymptotic  forms  indicated  in  Figure  12  follow  directly 
from  assuming  * - 0 as  |x  I — » and  writing 


Bound  states  of  the  scattering  problem  require  a(y')  = 0 and  Im(y')  > 0. 

The  corresponding  bound  state  eigenvalues  represent  the  kinks  and  breathers 
to  be  found  in  the  evolution  of  4>(x,t).  It  is  easily  shown  that  under  (B-4) 
the  upper  (lower)  half  of  the  y plane  maps  to  the  upper  (lower)  half  of  the 
y' -plane;  thus  the  condition  for  a bound  state  is  the  same  in  laboratory  (x,t) 
coordinates  as  in  the  characteristic  (£,t)  coordinates  usually  used  to  dis- 
cuss the  ISTM  for  the  sine-Gordon  equation. 

Since  the  bound  state  eigenvalues  are  independent  of  time,  they 
can  be  determined  from  the  initial  conditions  (B-3)  at  t = 0. 
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In  Regional) 


1 . 


V x 


4*j  = e 


*2  = ° 


In  Region  (2)  (Proceeding  as  in  [ 17] ) 

ill 


(B- 5a, b) 


4 = 


where 


(cos  mx  - sin  mx) 
2m 


- — sin  mx 
4m 


4m 


sin  mx 


(cos  mx  + rf-  sin  mx) 


10 


20 


(B-6) 


(2m)2  = (v')Z  + (\)Z 


(B-7) 


The  boundary  conditions  at  x = -p  determine  the  constants  and 


*20  93 


i , 

ly'  2 V’p 

4/,.  =(cos  mp  - sin  mp)e 
10  m 


■A  »,• 


*20  ~ 4m  Sln  mP)e2 


VP 


(B-8a,  b) 


Then  the  boundary  conditions  at  x = -fp  determine 

^ V *P  2V22 

a(y',  0)  = e [ (cos  mp  - sin  mp)  - (—  ) sin  mp]  (B-9) 


4m 


as  the  incident  wave  amplitude  at  time  t = 0. 
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,~A 


The  bound  states  are  determined  by  the  condition  a(y',  0)  = 0 
which  can  be  expressed  as 


cot  2mp 


_ ixl 


(B  - 10 ) 


This  condition  together  with  (B-7)  is  just  the  pair  of  equations  studied 

previously  in  connection  with  the  problem  of  fluxon  propagation  [ 17]  . The 

only  difference  is  that  y has  been  replaced  by  y'  which  is  related  to  y 

by  (B-4).  Thus  we  can  immediately  state  that  all  the  roots  of  (B-10)  which 

lie  in  the  upper  half  of  the  y' -plane  lie  on  the  imaginary  axis,  and  the 

construction  of  Figure  11  in  reference  [ 17]  will  find  them.  In  particular  we 

note  that  2m  p « rnr  for  most  of  the  roots,  so  from  (B-7) 
n 


^n  * p v (¥)2  ‘ (nw)2 


(B  - 11) 


To  see  which  of  these  roots  correspond  to  breathers  we  must  return  to 
the  y -plane  and  use  the  results  of  the  preceding  appendix.  First  we 


write  (B-4)  in  the  form 


2 jv'  ±n/(v’)2+  1 J . 


(B-12) 


As  V is  increased  from  zero,  threshold  levels  (V  ) are  reached  at 

’ n 

which  new  zeros  of  a(y',  0)  appear  in  the  upper  half  of  the  y'-plane. 
Each  enters  at  the  origin  and  moves  up  the  imaginary  axis  of  the  y'-plane. 
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From  (B-12)  a corresponding  pair  of  zeros  enters  the  y -plane  at  y = ±1/2 
and  move  up  the  circle  Iy  I = 1/2  indicated  in  Figure  U.  Thus 
stationary  breathers  are  generated  Just  above  threshold  by  the  initial 
conditions  (B-3).  From  the  construction  in  Figure  11  of  [17]  , it  is  easily 
seen  that  the  threshold  condition  for  the  n**1  breather  pair  is 

Vn  = (2n-l)w/p  . (B-13) 

From  (B-ll)  we  note  that  as  p -•■  <»  most  of  the  zeros  lie  at 
Y'  = iV/2.  For  V < 2,  (B-12)  indicates  that  they  lie  on  the  circle 
Iy  I = 1/2  in  Figure  11  and  cluster  at 

P = sin‘1(V/2)  . (B- 14) 

As  p — 00  and  V<2  in  Figure  12,  we  have  $ ■*  0 every- 

where whereupon  (1.1)  reduces  to  the  pendulum  equation 

4>tt  = sin  4>  . (B-15) 

This  extended  oscillation  maybe  viewed  as  "gas"  of  a large  number  of 
synchronized  or  coherent  breathers. 

It  is  interesting  to  compare  the  energy  of  this  coherent  state  with 
the  total  input  energy  UQ.  From  the  Hamiltonian  density  (A-12)  and  the 
initial  conditions  (B-3),  this  energy  is 

UQ  = V2P  (B-16) 


which  can  be  portioned  between  i)  kinks  (of  which  there  are  none  when 
V < 2),  ii)  breathers  and  ill)  radiation.  Since  these  three  energy  components 
are  positive  definite  [ 38  ] , the  energy  of  the  N breathers 


U 


NB  1 V 


(B-17, 


An  upper  bound  on  U is  obtained  by  noting  that  the  approximation  (B-ll) 

NB 

underestimates  |y^  | so 

and  from  (B -14 ) the  energy  of  the  n^1  breather 

unb>8v7^W~-  ,b-wi 

From  (B-13)  the  number  of  breathers  (N)  is  the  largest  integer  less  than 

Vpht 

2tt 

Thus  there  are  no  breathers,  and  no  breather  energy,  when 

Vp<  ir. 


In  general  the  total  breather  energy 

UNB  ” L «vy,-  (5BI)2  . (B-201 

n=l 
For 


V p » 2if 

the  right  hand  side  of  (B- 20)  can  be  approximated  by  the  integral 

J V2p  ^ Jl-yZ  dy  = V2p  ‘ 

Thus  (B-20)  and  (B-17)  together  imply 


1 

I 

5 
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as  p — oo  with  V < 2.  For  a number  of  breathers  large  compared  with 
unity,  essentially  all  the  input  energy  goes  into  the  breathers.  If 
Vp  = 2ttN,  (B- 20)  implies  uNg/u0  is  greater  than  55%  for  N = 2, 
greater  than  71%  for  N = 3 and  greater  than  79%  for  N = 4. 

Finally  let  us  suppose  that 

V » 2 

and  we  are  interested  in  the  behavior  of  the  dissipative  sine-Gordon 
equation  (1.2).  The  time  average  of  sin  <t  will  be  small,  and,  for 
sufficiently  large  p,  the  x derivatives  can  be  neglected.  Then  (1.  2) 
reduces  to 

dV 

“ * -«V  . (B-21) 

From  (B-13),  a breather  forms  every  time  V decreases  by  2ir/p. 

Thus  a rough  estimate  for  the  number  of  breathers  formed  per  unit 
time  and  per  unit  distance  is  given  by 


APPENDIX  C:  DOUBLE  PHASE  SOLUTIONS 


This  appendix  is  a summary  account  of  the  theory  of  multiphase 
solutions  for  nonlinear  wave  equations  which  has  been  developed  over 
the  past  few  years  by  Ablowitz  [9  - 12]  . 

Suppose  we  have  a nonlinear  wave  equation  for  <j>(x , t) 

N(<t,)  = 0 (C-l) 

which  can  be  obtained  from  an  Euler  variation  of  a Lagrangian  density. 
Thus 


N(<t»)  = H L(*t,  -4,x,  <j>) 


(C-2) 


Suppose  further  that  <fc(x,  t)  can  be  written  as  a double  phase  function 


<j>(x,  t)  = $(0^  02) 


(C-3) 


such  that  L is  a doubly  periodic  function  of  0^  : nd  0^.  Then  choose 


61  = Wlt  " ^1X 


°z-^-kzx 


(C-4a,b) 


and  average  L over  both  periods  to 

2tt  2tt 
' / 

4tt“  0 0 * "1  “ "2 


1 ’~2  / / w“lVVe  - -*l6Ve; 


(C-5) 


= H<v  u>2’  kr  k2)  ’ 
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If  w2’  ki  anc*  ^2  are  n0t  constant  but  allowed  to  vary  slowly 


with  x and  t,  two  dynamic  equations  for  this  slow  evolution  can 
be  obtained  from  Euler  variations  of  r with  0^  and  0^.  Thus 


- = 0 =>  = 0 


60 


'awj't  akj'x 


77T  £ = 0 =>  (T2*- L - = 0. 


60 


8“2t  9k2X 


(C-6a) 

(C-6b) 


Two  additional  equations  are  conservation  of  periods  for  the  two 

* 

components  or 


<“l>x  ♦ (kj)t  « 0 


(t°2)x  + «k2,t  = ° ’ 


(C-6c,d) 


A necessary  condition  for  this  description  to  be  valid  is  that  the  original 
equation  (C-l)  have  the  double  phase  solutions  indicated  in  (C-3)  for 
which  L is  doubly  periodic.  It  has  recently  been  demonstrated  that 
multiple  phase  solutions  [<f>(x,  t)  = #(0j,  ©2,  . . . , 0^)  for  any  finite 
integer  N]  exist  for  those  particular  nonlinear  wave  equations  which 
display  soliton  behavior  [13-15]  . This  class  includes  the  sine-Gordon 
equation  (1.1). 

A Lagranglan  density  for  (1.1)  is 

L = 2 4>x  ' 2 ^t  ' cos  <t>  • (C-7) 
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aim 


If  <t>(x,  t)  has  the  doubly  periodic  form  in  (C-3) 


L = l (k,*e  <■  k2%  )2  - \ * "2*0  ' ' cos  * 

1 z 1 z 


(C-8) 


so  (C-6a,b)  become 
2ir  2w 

A 
at 


2it  2tt 


/ /<“i*e1+"2*e1*62ldelci62ti;  / / <ki*e1tk2*e1%2,deide2=0 

2ir  2tr  2w  2v 

-2-  f f (w  +w,®  * )de,de^  + -*-[[  (k  +k.4>  ® 

9t  Vo  2 62  1 G1  02  1 2 8xJ0  0 2 62  1 9!  0 


)de1de2  = 0. 


(C-9a,b) 


These  are  (3.  22)  of  reference  [9] . 

Now  consider  (1.  2)  in  which  the  sine-Gordon  equation  is  made 
dissipative  through  addition  of  the  small  loss  term  Following  the 

discussion  in  Whitham  [7,  p.  510) , we  note  that  (C-9a)  can  be 
interpreted  as 


<*tlt  • <*x>lx  + <sln*>l  * ° 

where 


(C-10) 


2n  2ir 

>1  * / / < ■ 1* 


eideide2 


(C-ll) 


Noting  the  correspondence  between  (C-10)  and  (1.1)  in  the  form 


- <b  + sin  <t>  = 0.  we  can  perform  a similar  averaging  of  (1.  2) 
tt  xx 

to  obtain 


<+t>lt-<+x>lx+<sln  = 


(C-12) 


Similar  treatment  with  respect  to  the  averaging 

Zir  Ztt 

<‘>2zJ  f <-)»B  d8,de- 

2 0 0 ®2 

leads  to  the  pair  of  dynamic  equations 

2n  2ir  _ 2^  2ir 


(C-13) 


«/  / (“i®o,+"2*e,*9,|deide2  + ^/  / <ki*6.+k*9  *e Jdelde2  • 


0 0 


'1  - V1  "2 


0 0 
2lT  2lT 


1 1 2 


= + -2®e  • e l<*V° 


0 0 


1 2 


= -»  f V VeI+“2*oI*eJ)deldez  ' 

(C-I4a,b) 

These,  together  with  the  conservation  equations  (C-6c,d)  are  the  dynamic 
equations  used  in  this  paper.  In  Section  III,  only  a single  phase  is 
assumed  (3.3)  so  (C-I4a)  reduces  to  (3.  5a).  In  Section  IV,  both  phases 
are  employed  for  the  d.c.  components.  It  is  necessary  in  either  case  to 
assume 
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a « 1 


so  that  the  dissipative  variation  over  a cycle  is  small  and  the  averaging 
indicated  in  (C-5)  is  well  defined. 

Ablowitz  [9]  has  suggested  a way  to  get  some  integral  constraints 
on  the  function  ®(6j,  0^).  Considering  (C-3)  as  a change  of  independent 
variables,  (1.1)  becomes 


Veje/  2h®e1e2t92*e2e2+sin®  =0 


(C- 15) 


where  g^  s and  h = Wj«2  - e<luat^on  can  be  written 

in  the  form  of  a conservation  law  as 


9l*e1  - i 92®e2  - cos®>  + ^ <h*el  4 = 0 


or  as 


(C-16a,b) 


— ^ — (h  ^ + g 1 ® ) + — ^ — (~  g 4>  — - g,®  ^ - cos  ® ) = 0. 

8G1  02  91  G1  Q 2 dQz  2 92  G2  2 91  Gj  ' 


From  these  it  follows  that 
2 ir 


/ <lVo1-I92®e2-cos®we2-Ei 


and 


2" 


(C-I7a,b) 


/ (2  92®62  " 2 *1®^  ' °°S  ®ldei  = "E2  ' 


If  ® is  independent  of  62,  (C-I7a)  is  satisfied  by  the  elliptic  function 
in  (3.4). 
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Under  the  assumption  h = 0,  (C-15)  becomes 


g,<X>  + _ + sin  9 = 0 

1 Vl  2 e2e2 


which  Lamb  [ 36]  has  shown  to  have  solutions  of  the  form 


* = 4 tan  [ f(0j)g(02)] 


(C-19) 


where  f and  g are  elliptic  functions.  This  is  the  doubly  periodic 
function  suggested  recently  by  Ben-Abraham  [ 42]  . For  our  purposes, 
however,  the  assumption  h = 0 is  too  severe.  It  requires 


f CO  CO  i 

'_1\  -2  _ 1 

k k _ 1 

II  2/ 


(C-  20) 


which  means  that  the  two  component  traveling  waves  must  have  phase 

velocities  in  the  same  direction  with  one  going  faster  than  the  characteristic 

velocity  of  (1.1)  and  the  other  slower. 

From  references  [13-15]  it  seems  that  a simple  prescription  for 

constructing  a more  general  double  phase  solution  to  (1.1)  would  be  as 

follows:  i)  Express  the  elliptic  function  defined  in  (3.4)  in  terms  of 

theta  functions,  ii)  Augment  the  arguments  of  these  theta  functions  from 

0 to  6 + 0 , and  iii)  Check  that  the  resulting  expression  is  a solution 
1 2 

of  (1.1). 
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